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Abstract 



We propose a novel approach for studying i/^ — > Uj- oscillations with ex- 
tragalactic neutrinos. Active Galactic Nuclei and Gamma Ray Bursts are 
D ■ believed to be sources of ultrahigh energy muon neutrinos. With distances of 

■ 100 Mpc or more, they provide an unusually long baseline for possible detec- 

, tion of Ufj, — > Vr with mixing parameters Am^ down to 10"^'^ eV^, many orders 

of magnitude below the current accelerator experiments. By solving the cou- 
pled transport equations, we show that high-energy i^t-'s, as they propagate 
through the earth, cascade down in energy, producing the enhancement of the 



0^ ■ incoming u-j- flux in the low energy region, in contrast to the high-energy f^'s, 

which get absorbed. For an AGN quasar model we find the Vr flux to be a 
I factor of 2 to 2.5 larger than the incoming flux in the energy range between 

' 10^ GeV and 10^ GeV, while for a GRB fireball model, the enhancement is 

10%-27% in the same energy range and for zero nadir angle. This enhance- 
ment decreases with larger nadir angle, thus providing a novel way to search 
^ . for z^T- appearance by measuring the angular dependence of the muons. To 

D i illustrate how the cascade effect and the Uj- final flux depend on the steepness 

' of the incoming z^,-, we show the energy and angular distributions for several 

generic cases of the incoming tau neutrino flux, Fj^ ~ for n = 1,2 and 
^ ■ 3.6. We show that for the incoming flux that is not too steep, the signal for 

^ . the appearance of high-energy i^r is the enhanced production of lower energy 

fi and their distinctive angular dependence, due to the contribution from the 
r decay into fi just below the detector. 

Recent Super-Kamiokande (SuperK) measurements of the low atmospheric z/^/z/g ratio 
and the strong zenith angle dependence of the z/^ events [0 suggest oscillations of z/^ into 
z/^ with the parameters sin^ 26' > 0.7 and 1.5 x 10~^ < Am^ < 1.5 x 10"^ eV^ This is 
in agreement with previously reported results on the atmospheric anomaly by Kamiokande 
0] and MACRO and is consistent with limits from other experiments, e.g., CHOOZ [Q. 
Confirmation of — > Ur oscillations and determination of neutrino mixing angles would 
be a crucial indication of the nature of physics beyond the Standard Model. The firmest 
confirmation of this hypothesis would be via detection of r leptons produced by charged 
current interactions of z/^'s resulting from oscillations of z/^'s, which is extremely difficult 
with current neutrino experiments. 

In this letter, we propose a study of z/^ Vr oscillations with extragalactic neutrinos. 
Large volume neutrino detectors and the prospect of astrophysical neutrino sources put Vr 
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detection in the realm of possibility. The large distances involved for astrophysical sources, 
on the order of one to thousands of Megaparsecs, make the next generation of neutrino 
experiments potentially sensitive to neutrino mass differences as low as Am^ ~ 10"^"^ eV^ 
IP). Over such long baselines, half of the neutrinos arriving at the earth would be z/^'s in 
oscillation scenarios, the other half being z/^'s. By observing both u^^ and from extragalatic 
sources such as Gamma Ray Bursts (GRBs) and Active Galactic Nuclei (AGN) [0, 
neutrino oscillation hypothesis would be confirmed and models of these sources would be 
tested. 

The effect of attenuation of the neutrino flux due to interactions of neutrinos in the 
Earth is qualitatively different for z/^ and u^-. Muon neutrinos are absorbed by charged 
current interactions, while tau neutrinos are regenerated by tau decays. The Earth never 
becomes opaque to Ur, though the effect of z/^ — ^ t — * z^r interaction and decay processes 
is to degrade the energy of the incident Ut-. The identical spectra of z/^ and u^- incident on 
the Earth emerge after passage through the Earth with distinctly different spectra. The 
preferential penetration of z/,- through the Earth is of great importance for high energy 
neutrino telescopes such as AMANDA, NESTOR and ANTARES. 

We consider z/^ and Ur propagation through the Earth using a similar procedure to the 
one outlined for z/^ in Ref. [Q. We show that the energy spectrum of the becomes enhanced 
at low energy, providing a distinctive signature for its detection. The degree of enhancement 
depends on the initial neutrino flux. We consider initial fluxes F° ~ E^'"' for n = 1,2, 3.6, a 
GRB flux and an AGN flux [0 . We solve the coupled transport equations for lepton and 
neutrino fluxes as indicated below. 

Let Fy^{E,X) and Fr{E,X) be the differential energy spectrum of tau neutrinos and 
tau respectively at a column depth X in the medium defined by 



where p{L) is the density of the medium at a distance L from the boundary measured along 
the neutrino beam path. Then, one can derive the following cascade equation for neutrinos 

as. 
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and for taus as, 
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dFriE,X) ^ F,{E,X) _ F,{E,X) 
dX K{E) pf-{E,X,e) 

^{u^N^TX;Ey,E). (2) 

The first term in Eq. (1) is a loss due to the neutrino interactions, the second is the 
regeneration term due to the neutral current, the third term is a contribution due to the tau 
decay and the last term is the contribution due to tau interactions. 

In Eq. (2), the first term is a loss due to tau interactions, the second term is a loss 
due to the tau decay, while the last term is a contribution from neutrino charged current 
interactions. As a practical matter, tau decays are more important than tau interactions at 
the energies considered here, though interactions become more important at higher energies 
0. We neglect the tau interaction terms in Eqs. (1) and (2) in what follows. 

Here X{E) is the interaction length and pf^'^{E,X,9) is the decay length for tau. They 
are defined as. 
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where Nt is the number of scatterers T in 1 g of the medium, o"*^(-E) is the total cross 
section for the uT interactions and the sum is over all scatterer types {T = N,e, . . .), (r is 
the mean lifetime of tau and g is the density of matter in the earth. Scatterings of neutrinos 
and taus with nucleons (A^) are most important, so we approximate 



where Nq is Avogadro's number. We use the CTEQ5 parton distribution functions to evalu- 
ate neutrino cross sections [0. We have previously calculated charged and neutral current 
energy distributions, dn/dE, and the total cross section, a^"^{E) taking into account recent 
improvements in our knowledge of the small-x behavior of the structure functions |TT]. 

To simplify the solution to the equation for the tau flux, we approximate the X and 
9 dependent density of the earth by the average of the density along the column depth of 
angle 9: 

let us define the effective absorption length A^{E,X) by 

X 



Following Ref. 



F,{E,X) = F:^{E)exp 



ME,X) 



(3) 
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It is convenient to define 



K{E,X) 



1-Z,{E,X) 



(4) 



where Z,y{E,X) is a positive function (we will call it Z factor in analogy with the hadronic 
cascade theory) which contains the complete information about neutrino interaction and 
regeneration in matter. 

Assuming that there is no significant contribution to the neutrino flux from decaying 
particles (as would be the case for muon neutrinos) using the above equation we can find an 
implicit equation for Z from the transport equation p|, 



where 
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and da^N^uxiy, E)/dy is the differential cross section for the inclusive reaction uN vX 
(with Ey the incoming neutrino energy and y the fraction of energy lost) and Ey = E/{l — y). 

Naumov and Perrone have shown that by iteratively evaluating Eq. (5), starting with 
Z^^^ = 0, the solution for muon neutrinos quickly converges for a wide range of starting 
fluxes. 

By a similar procedure, the coupled differential equations for tau neutrinos including tau 
production and decay can be iteratively solved. The tau flux generated by charged current 
interactions including the loss term due to its decay is 
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The Z factor for the tau neutrino flux is then 

Z = Zi, + Zj- 

where Z^, is given by Eq. (5) and 
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We include decay modes in ^ff'^{y, E) as in Ref. [T^ and a constant energy distribution 
for the remaining branching fraction not included there. In Eqs. (5) and (8), the Z factors 
implicit in Aj, are Z = Z^y + Zt-. 

In the iterative solution of the equation for Z, one has the option of picking the initial 
value Z^^\ We have chosen the X and E dependent solution to the cascade equation for 
the z/^ flux, namely the solution to Eq. (5). 

To demonstrate the importance of regeneration of tau neutrinos from tau decays, we eval- 
uate the tau neutrino flux for several input neutrino spectra and compare to the attenuated 
z/^ flux. For the incoming neutrino spectrum we use 

F^AE)=K(^yj^), (9) 



where K, n, Eq and Ecnt are parameters and is a function equal to at t > 1 and 1 at 
t <Cl. We use (p{t) = 1/ [1 + tan (7rt/2)] (t < 1) and E^nt = 3 x 10^° GeV and = 1 PeV. 
For n = 1, we introduce a smooth cutoff by multiplying Eq. (9) by a factor {1 + Eq/E)^'^ 
and by setting Eq = 100 PeV. To evaluate the depth as a function of nadir angle, the density 
profile of the earth described in Ref. |11| is used. 



In Fig. 1 we show the nadir angle dependence of the ratios of the fluxes calculated via 
Eqs. (1-8) to the input flux F^{E). All fluxes are evaluated as a function of nadir angle, at 
the X value for the surface of the earth and for energies 10^ GeV, 10^ GeV and 10® GeV for 

and z/t- assuming F^{E) given by Eq. (9). For an incoming flux, n = 1, we find that the 
z/t- flux is enhanced (relative to the incoming v-r flux) for all nadir angles for E^^ = 10'^ GeV 
and 10^ GeV, and for 9 > 20° for E^^ = 10® GeV. The peak of the enhancement gets shifted 
toward the higher nadir angles as the energy increases. This is due to the fact that high 
energy z/,- can remain high energy if the column depth is small, i.e. for large nadir angles. 
In case of the steeper incoming flux, n = 2, we find that z/^'s are less attenuated than the 
z/^'s, and the expected enhancement at low energy is not evident due to the steepness of the 
flux. For even steeper flux, n = 3.6, the difference between z/,- and z/^ flux is very small. 

In Fig. 2 we show the energy dependence of the the ratio of fluxes for nadir angles ^ = 0, 
6 = 30° and 6 = 60° for z/^ and Ur with F^{E) given by Eq. (10). For small nadir angles, 
6 = and 30° and F^{E) ~ 1/E we find that enhancement of tau neutrinos is in the energy 
range of 10^ GeV and 10^ GeV, while for 6 = 60°, the enhancement extends up to 10® GeV. 
In contrast the z/^ flux is attenuated for all the nadir angles. When the incoming flux is 
steeper, n = 2, the Ur flux appears to be attenuated at high energies, although less than 
the z/^ flux. For n = 3.6, the energy dependence of these two fluxes is very similar, they are 
both reduced at high energies, and the effect is stronger for smaller nadir angle, since in this 
case the column depth is larger and there are more charged current interactions possible. 

In case of the AGN quasar model 0, we find that the z/^- flux is a factor of 2 to 2.5 times 
larger than the input flux, for nadir angle, 6 = 0. This is shown in Fig. 3. For larger angles, 
the effect is smaller. Detection of AGN neutrinos would be optimal for small nadir angles 
and for z/,- with energy of 10^ GeV to 10^ GeV. 

We also present results for the u^- flux for the case of GRB fireball model 0. We find 
that due to the steepness of the input flux for E^,^ > 100 TeV, the u^- flux is enhanced only 
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by about 10 — 27%, depending on the energy and nadir angle. We show the energy spectrum 
of the ratio of z/^ flux to the input flux F° in Fig. 4. 

We expect that the next generation of neutrino telescopes will be able to detect the high 
energy neutrinos from AGN and GRBs. We have previously shown that most of the ex- 
tragalactic neutrino fluxes exceed the atmospheric neutrino background for neutrino energy 



greater than ~ 10 TeV which may enable the detection of the extragalactic neutrinos |]TT 
A search for high-energy z/^ appearance would look for enhanced muon rates. 

The enhanced rates of muons come from muonic decays of the r produced by charged 
current interactions in or near the instrumented detector volume. The muon rates would 
be enhanced at low energy (~ 10 — 100 TeV) and for small nadir angles. In the case of 
an AGN quasar model, we find that the v^- flux enhancement is very distinct, while for a 
GRB fireball model, the effect is only 10 — 27%. This is due to the fact that GRB input v-r 
flux is much steeper, thus there are not many high-energy i/^-'s that would contribute to the 
enhancement at low energy. We have also shown that the low energy pile up is significant 
only for incoming fluxes that are less steep than in the high energy region. For an 

incoming flux which is proportional to 1/-E, we find the angular distribution of v^s to be 
significantly different than for the f^'s. 

We have proposed a novel way of detecting appearance of extragalactic high-energy Vt 
by measuring the angular and energy distribution of muons with energy above 10 TeV. This 
would give an experimental signature of z/„ oscillations with Am^ as low as 10~^^ eV^ 
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FIGURES 
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FIG. 1. The nadir angle dependence of the ratio of fluxes for energies 10^ GeV, 10^ GeV 
and 10^ GeV for the z/^ and z/^ assuming F^{E) ~ E''" with n = 1,2 and 3.6. 
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FIG. 2. The energy dependence of the ratio of fluxes for nadir angles 6 = 0, 6 = 30 and 
6 = 90° for z/^ and z/,- assuming F^{E) ~ with n = 1,2 and 3.6. 
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